Superconductivity in multi-orbital t — Ji — J2 model and its impli- 
cations for iron pnictides 
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Abstract. - Motivated by the bad metal behavior of the iron pnictides, we study a multi-orbital 
t — Ji — J2 model and investigate possible singlet superconducting pairings. Magnetic frustration 
by itself leads to a large degeneracy in the pairing states. The kinetic energy breaks this into 
a quasi-degeneracy among a reduced set of pairing states. For small electron and hole Fermi 
pockets, an Aig state dominates over the phase diagram but a Big state has close-by energy. In 
addition to the nodeless Aig s^2y2 channel, the nodal Aig s^ij^yi and B\g d^2_y2 channels are 
also competitive in the magnetically frustrated Ji ~ J2 parameter regime. An Aig -\- iBig state, 
which breaks time-reversal symmetry, occurs at low temperatures in part of the phase diagram. 
Implications for the experiments in the iron pnictides are discussed. 



Introduction. — The discovery of high temperature 
superconductivity in the iron pnictides [lH3| has spurred 
tremendous experimental and theoretical interest in these 
systems. While they are not Mott insulators, the un- 
doped iron pnictides are nonetheless "bad metals". This 
fact has motivated the placement of these systems in an 
intermediate coupling regime close to the boundary be- 
tween Mott localization and itinerancy [SHZ|, where the 
Coulomb interactions bring about non-perturbative effects 
in the form of incipient lower and upper Hubbard bands. 
Accordingly, the low-energy Hamiltonian contains quasi- 
localized moments with Ji — J2 superexchangc interac- 
tions [6HT2] , and it corresponds to an effective multi-band 
t — Ji — J2 model for the carrier-doped systems. Impor- 
tantly, Ji - J2 > Ji/2 PlIMTUlfn] so that the system 
not only has a (tt, 0) antiferromagnetic order, as seen ex- 
perimentally [4], but also exhibits strong magnetic frus- 
tration. The incipient Mott picture is supported by the 
observations of the Drude- weight suppression p^HIB] . as 
well as the temperature-induced spectral-weight transfer 
[151 [171 [18] . In addition, the quasi- localized moments are 
supported by the inelastic neutron scattering experiments; 
Ref. [H] , for instance, observed zone-boundary spin waves 
and showed that Ji is indeed comparable to J2. (By 
contrast, Fermi surface in the magnetically ordered state 
does not directly probe the strength of electron correla- 
tions [7].) Alternatively, perturbative treatments [201121] 



of the Coulomb interactions have been used to study the 
magnetism of the iron pnictides. The weak-coupling ap- 
proaches, mostly based on spin fluctuations, have been 
extensively used to address the superconductivity [^IH^B] . 
By contrast, strong-coupling studies of superconductivity 
have been more limited [27l - l30] . 

Experimentally, the pairing symmetry in the iron pnic- 
tides has remained inconclusive. The angle resolved 
photoemission [311121] (ARPES) and the Andreev spec- 
troscopy [33] results suggest a nodeless gap. In contrast 
the nuclear magnetic resonance [34] (NMR) and some pen- 
etration depth |35| measurements suggest a nodal gap. 
The experimental results seem to vary among pnictidc 
compounds. As an example, the P-doped BaFe2As2 ap- 
pears to have nodal gaps [3S1I37] and this is in strong con- 
trast with its K-doped counterpart. 

In this Letter we address the pairing in the pnictides 
from the incipient Mott approach, and treat the competing 
pairing channels on an equal footing. We are particularly 
motivated to consider the effect of the Ji ~ J2 magnetic 
frustration, and show that it leads to a quasi-degeneracy 
among several paring states. 

The Hamiltonian is given by 



H = Ht+H 
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H 



J2- 



(1) 



The kinetic part is Ht = ~ ^ 



i<j,a,l3,s ''ij ^ias'^jfi 



h.c 



where 



creates an electron at 
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site i, with orbital a and spin projection s; ^ is the 
chemical potential and t"^-^ the hopping matrix. The 
nearest- neighbor (n.n., (ij)) and next-nearest-neighbor 
(n.n.n., {{ij})) exchange interactions are, respectively, 

Hji ~ J2{i]),a,p i^^ia • Sj(3 - \nianjp^, and Hj2 = 

kT.s,s' 4as^ss'C^as' and n,^ = 'Es4as(^^c,s, with CT rep- 
resenting the Pauli matrices operating on the spin indices. 
The above Hamiltonian is augmented by the appropriate 
occupancy constraint for the fermions. Note that, while 
the nearest-neighbor interaction has been evidenced to be 
spatially anisotropic in the (tt, 0) coUinear antiferromag- 
netic state [9l[T9], it is expected to be isotropic in the 
tetragonal paramagnetic phases and this is consistent with 
spin dynamical measurements [38(139). 

Two-orbital model. — We will first consider a two 
orbital model, retaining only the dxz{c^ = 1) and dyz{a = 
2) orbitals, and later we consider a five orbital model to 
better address the fermiology [2Tl[26l|4^. We further as- 
sume J"'^ = Ji5a,p, and > 0, and will subsequently 
address the role of interorbital exchange couplings. 

Under the tetrahedral point group symmetry trans- 
formations (D^fi), dxz and dyz orbitals transform re- 
spectively as X and y coordinates. The kinetic en- 
ergy part of the Hamiltonian is invariant under all 
point group symmetry operations {Aig) and has the 
following form in the extended Brillouin zone, Hq = 
I]k,s V'L Kk+To + Ck-Tz -I- (,\^xyTx\ ipUs, whcrc k = {kx, ky) 

and tjjks — ("^kis' '^k2s)- "^^^ identity and Pauli matri- 
ces (tq, Ti) operate on the orbital indices, and = 

— {ti + t2){cos kx + cos ky) — At^ cos kx cos ky — fi, ^k- ~ 

— {ti — t2){coskx ~ cos ky), ^kz:y = — 4^4 siu fcj, siu fcy are 
respectively Aig, Big, B2g functions. The band disper- 
sion relations £k± = Ck+ ± \l^k- + ^L-y §i'^^ ^'^^^ 
electron pockets at k = (7i',0) and (0, tt), and two hole 
pockets at k = (0,0) and (7r,7r). The carrier doping 
(5=1^^ nia—2\. Mostly we use the minimal tight-binding 
model of Ref. gT], ti = -t, t2 = 1.3<, = t4 = -0.85t 
obtained from a fitting of the LDA bands. 

Magnetic frustration and degeneracy of pairing states. 
We start from the case with a vanishing kinetic energy, 
in order to highlight the connection between magnetic 
frustration and enhanced degeneracy of pairing states. 
We define the intra-orbital spin-singlet pairing opera- 
tors Ae,Qa = (CiatC-t+eaJ, - CialCi+eat) f^, whcrC G = 

x,y,x i y. Without the kinetic term, the problem de- 
couples in the orbital basis and we can drop the or- 
bital indices. When Ji dominates, two degenerate pair- 
ing states Sx2+y2 and dx^-y^, respectively defined by the 
pairing functions 53,2 ±y2 k = cosfc^, ± cos ky are naturally 
favored. In real space, they respectively correspond to 
Ax = ±Ay = Ao, with Ax+y = Ax-y = 0. These 53,2+^2 
and dx^-y^ states are degenerate because the symmetry 
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Fig. 1: Zero temperature phase diagrams of (a) a two orbital 
model and (b) a five-orbital model, both for electron doping 5 = 
0.14. The onset of Big and B2g phases are respectively marked 
by black and red solid lines. The dotted line characterizes a 
cross-over between s 2^2 [region (1)1 and s 2", 2 [region (2)1 as 

X y X -\-y 

the dominant component in the Aig pairing. 

operation Cmx+ny ^^'"^"^^^^^^'^Cmx+ny transforms them 
into each other. 

When J2 dominates, the Sx2y2 and dxy states, re- 
spectively defined by the pairing functions gx2y2^]^ = 
cos{kx + ky) + cos(kx — ky) and gxyM = cos{kx — ky) — 
cos{kx + ky), are preferred, and they are degenerate. In 
real space, they correspond to Ax+y = iAx-y = Aq, 
with Ax = Ay = 0. The S2;2j^2 and dxy states trans- 
form into each other by the following symmetry opera- 
tion: We break the square lattice into two interpenetrat- 
ing sublattices;on the even sublattice (m + n = even), 
Cmx+ny ^ e*("+"+i)''/2cm£+„ji , and on the odd sublat- 

tifP r - , - —5. pi(in-n+l)Tr/2 . . 
Lii^c, ^mx+ny ~ ^ ^mx+ny • 

As we tune the ratio J2/ Ji, we expect a level crossing in 
the magnetically frustrated regime, J2 ^ Ji- We can then 
anticipate that magnetic frustration promotes an enlarged 
degeneracy among the Sx2y2, Sx2+y2, dxy, and dj.2_y2 pair- 
ing states. 

The effect of the kinetic energy. When the ki- 
netic term is incorporated, it lifts the exact degenera- 
cies of the paired states discussed above. We study 
the full problem using a mean-field decoupling |42) of 
the two-band t — Ji — J2 model. To set the stage, 
we note that the D^h point group symmetry opera- 
tions allow the following four classes of pairing states 
for an orbitally diagonal Ji — J2 model [53]: (i) Aig : 

[S^2^y2gx2 + y2M + Sx2y2gx^y^M'^0 + "j.2 _ 5a;2 -y2 ^k^z ; (,nj 

B7-Blg . r Big Big 1 

Ig ■ "^2_y2 32:2— y2 kTo -T [S ^2 _^y2 gx^ +y^ ,k ~^ S ^2y2 gx^y^ ,l(.\'''z] 

A B 

(iii) A2g : dxy" gxy,kTz; and (iv) B2g : dxy'' gxy^TQ. Each 
pairing channel will have different symmetry depending on 
whether it is associated with tq or in the orbital space; 
this distinction is denoted by the superscripts. The eight 

Al 

pairing amplitudes s^2!^y2 etc. arc linear combinations of 
eight intra-orbital pairing amplitudes Ae.aa. 

Mean field theory. We decouple the exchange interac- 
tions in the pairing channel. In terms of = (V'kf i ^-k^). 
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Fig. 2: (Color online) The amplitudes of different pairing gap 
components of a two orbital model for electron doping 5 — 



0.14. For J2 > Ji and J2 < Ji, 



and d 



respectively the dominant pairing channels. 



the Hamiltonian becomes 



mf 



At 



-hk 



(2) 



where hk = Ck+To + ^k-Tz + ^kxyTx, and Ak = 
diag[Ak,ii, Ak,22] is the orbitally diagonal gap matrix, and 
Ak,Qa = J2e Je^e,aa cos(k • e), with Je = Ji for e = x,y, 
and Je = J2 foi' e — X ± y. Diagonalizing Hmf we obtain 
the quasiparticle dispersion spectra i?k,±- We determine 
the pairing gap matrix by minimizing the ground state 
energy density 



(3) 



with respect to all Ae^aa- The primary effect of the con- 
straints is to renormalize the kinetic energy via (both or- 
bitally diagonal and off-diagonal) t — > tS/2, where t is 
the kinetic energy scale. Our results, with an implicit 
treatment of the constraint through a band renormaliza- 
tion, remain largely unchanged when the constraints arc 
explicitly incorporated. 

When kinetic energy is absent, the s.j;2^y2 and d^2_y2 
states are indeed degenerate, and each has a ground state 
energy ^ — 0.17Ji. Likewise, the energy of either 53,2^2 or 
dxy state is « — O.I7J2, and all four paired states become 
degenerate exactly at Ji = J2, as anticipated earlier. 

We now turn to the results for the full problem in the 
presence of the kinetic terms. Wc show an illustrative 
zero-temperature phase diagram for < Ji, J2 < 2i, in 
Fig. 1(a) corresponding to an electron doping S 
A 



0.14. 

J2 plane. In a sizable 



iig pairing exists in the entire Ji 
portion of the phase diagram, the pairing is in the pure 
Aig class. In this region s^s^^s is the dominant pairing 

channel and coexists with the subdominant d 2" 2 chan- 

^ y 

nel. The onset of Big pairing state is marked by a solid 
line, which corresponds to a second order phase transi- 



tion, d 



Big 



-V 



and 



are respectively the dominant and 




Fig. 3: (Color online) The pairing amplitudes for d^^, d^z or- 
bitals, obtained from a five orbital model for 5 = 0.14. For 
Ji'^ the dominant pairing channel is s^2^^2 for 

J 2 <^ Ji, d^2^ 2 s-nd s^l^ 2 pairing channels are dominant 
and nearly degenerate pairing channels. Compared to the two- 

A-i 

orbital model s^2^ 2 channel is more competitive and becomes 
significant for relatively smaller value of Ji. 



the subdominant components of Big phase. In the large 



Ji limit (to the right of the dotted line) , s 



2;2_j_y2 



becomes 



the dominant component of the Aig phase; the dotted line 
represents a crossover. The phase diagram for hole doping 
is identical, but the onset of Big pairing occurs for larger 
Ji. 

The competition among different symmetry classes, and 
the nature of the ground states are demonstrated in Fig. [51 
which plot the pairing amplitudes as a function of the 
coupling constants. The degeneracy between the 3^.2^,2 
and dxy channels, occurs at very large J2 limit and is 
not showed in Fig. [5] For moderate values of Ji , J2 , 
quasi-degeneracy among the pairing states occurs in the 
magnetically- frustrated region corresponding to J2 Ji- 
In this region, Fig[5]shows that the weights of the s^2'^2 and 

d^2^y2 components are comparable, and those of d'^^'Ly^ 

and s^2^^2 are close-by. This quasi-degeneracy also under- 



lies the phase diagram shown in Fig. 1(a) By contrast, for 



moderate values of J2 and with J2 S> Ji , s 2 "2 dominates 

_ y 

Big 

over s 2 2 



,A2g ,S2g 

„2 , iJ-xy , fJ-xy 



channels; likewise, for moderate val- 



]3 A 

ues of Ji and with Ji ^ d^2'L--^ dominates over d 

A-lg 
S 2 , 

x^+y 

At the T 



2_,,2 : 



2' ^x'^+y^ 



channels. 

limit we study, the co-existing pairing 
channels will lock into a definite phase. Consider first 
the case of coexisting s^i^^ and d^i'^ 2 ■ We define Ak = 

sf2y2ffa:2y2_kro + (if2!.y2 5a;2-2/2,kT"z, and use to denote the 
relative phase. Inserting this into Eq. ([3]), we find that 
(j) ~ TT corresponds to the ground-state energy minimum, 
with df /d(j) = and df /d(j) > (0 = is a ground-state 
energy maximum). The Aig gap is therefore real. More- 
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over, the relative minus sign associated with (j) = tt imphes 
that the pairing function, in the band (as opposed to the 



orbital) basis, is sJy2gx^y'',kTo - cL^'^Lyi \ 

^kxyTa;)(Ck- +^kxy)~^^^- Hence intrabaud pairing function 
changes sign between the hole and electron pockets near 
k — (0,0) and k = (tt, 0), (0, tt) respectively. A similar 
argument shows that the relative phase between d^2'!_y2 
and s^2!^n is tt, and the gap changes sign between the two 
electron pockets near k = (tt, 0), (0, tt) respectively. 

Consider next the case of coexisting s'^i^i and d^i" 2 ■ 
A similar analysis shows that, this time, (/> = 7r/2 corre- 
sponds to the ground-state energy minimum (while (j) — 
0,7r represent energy maxima). The resulting Aig + iBig 
phase describes the state to the right of the solid line in 
Fig. 1(a) This phase simultaneously breaks time rever- 



A19 



sal and four-fold rotational symmetries, but preserves the 
combination of the two symmetries. Such a state also oc- 
curs in a phenomenological Landau-Ginzburg theory j44j . 

The role of fermiology can be clearly illustrated in the 
linearized gap approximation. For a set of small Fermi 
pockets at (0,0), (7r,0) and (0, tt), compared to gxy,'k, 
gx^y^,k has larger overlap with the pairing kernel. Thus 
s^2y2 and s^2y2 gaps have higher T^s compared to dxy" 

and dxy" gaps; they become degenerate only in the large 
J2 limit. Similar reasoning shows that, unless a thresh- 

^2-y2 gaps have 
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old value for Ji is exceeded, d 



and d 
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higher T^s compared to the 53.2+^2 and 53.2+^2 gaps. Re- 
lated observations were made by Seo et al. [28]. If we 
consider a band structure that produces large pockets 
[27], dxy" replaces s^2^^2 as the dominant pairing state for 
J2 > Ji, and we observe a competition between d^y' and 
d^2^ 2 pairing states. Note that the magnetic oscillation 
and Irpes measurements suggest small sizes of the Fermi 
pockets for the iron pnictides, and this should make our 
Aig and Aig + iBig phases more feasible. 

Effects of inter- orbital exchange couplings. We stress 
that, while the detailed nature of the lattice symmetry and 
orbitals for the single-particle energy dispersion is impor- 
tant for a proper description of the Fermi surface, its coun- 
terpart for the exchange interactions is not obviously so. 
In the absence of a detailed knowledge about such struc- 
ture in the exchange interactions, we have considered the 
simplest description and focused on the associated prop- 
erties that arc qualitative and robust. For instance, we 
have so far considered intra-orbital exchange interactions. 
Inclusion of inter-orbital super-exchange interactions docs 
not change the phase diagram. Now the Aig phase will 
have a small inter-orbital dxy" component and the Big 
phase will remain unchanged. Hund's coupling may lead 
to on-site triplet pairing, but this is prevented by the on- 
site Coulomb repulsion that is built in our model. Inter- 
site triplet pairing may arise if the Hund's coupling is com- 
parable to Coulomb repulsion |29j . but this is unlikely on 



general or ab initio grounds; it is also unlikely on empir- 
ical grounds since experimental evidence has so far been 
overwhelming for singlet pairing. Using a multi-orbital 
Hubbard model as their starting points, both the strong 
coupling calculations of Ref. [55] and the RPA calculations 
of Ref. [U find that a moderate Jh enhances repulsive 
inter-electron-pocket pair scattering and leads to stronger 
Big pairing. Thus Jh can only enhance the Aig + iBig 
part of the phase diagram. We have preferred to consider 
the multiband t — Ji — J2 model, since Ji and J2 can 
be more readily connected with the magnetic frustration 
physics. 

Five-orbital model. — To understand the robust- 
ness of our two-band results, against the inclusion of ad- 
ditional bands, and better address the fermiology we have 
considered a five-band model with the kinetic terms ac- 
cording to Ref. [in]- To capture the important results 
within a simple model of interaction, we again choose 
a Ji — J2 model with J"^ = JiSap- Now the gen- 
eral intra-orbital pairing matrix has the form Ak = 
Eadiag[A« ii,A«22,A-33,A-44,A- 55], where the in- 
dex a corresponds to Sx2+y2, dx2-y2, Sx2y2 and dxy sym- 
metries. The results of minimizing the free energy with 
respect to all the twenty complex pairing amplitudes are 



given in an illustrative phase diagram Fig. 1(b) and in 



Fig. [S] which shows the competition among the pairing 
amplitudes for xz and yz orbitals. We find that the compe- 
tition between Aig and Big pairings is a robust effect. In 



contrast to the two band case d, 

„Alg 



Aig 



amplitude is reduced 



and s^2!ly2 amplitude is enhanced in the entire phase di- 
agram. Indeed, for the magnetically-frustrated Ji ^ J2 
region, the nodal Aig Sx2+y2 and Big dx2_y2 states are 
competitive against the nodeless Aig 3^.2^,2 state. Com- 
pared to the two band case the B2g pairing occurs for 
smaller J2, but it still occurs in the limit when J2 is big- 
ger than the kinetic energy scale. The other three orbitals 
also demonstrate similar competition among Sx2y2, 3^2 ^y2, 
dx2^y2, and d^y pairings. 

Experimental implications of quasi-degenerate 
pairing channels. — Some of the weak-coupling stud- 
ies [mim] have also indicated a competition between 
various pairing states. However, the weak-coupling ap- 
proaches are typically restricted to an instability analysis 
of the linearized gap equations. By contrast, the strong- 
coupling approach used here has the advantage of readily 
considering the non-linear gap equations. Our non-linear 
analysis is important in bringing out the connection be- 
tween the quasi-dcgeneracy in the pairing channel and the 
Ji J2 magnetic frustration. With only one require- 
ment (Ji ~ J2), which is linked to magnetic frustration 
and has been supported by both theoretical considerations 
[6ll8lI10|13] and inelastic magnetic experiments [T9ll38] , our 
result provides a parameter-insensitive mechanism for the 
near degeneracies of various paring states. Moreover, our 
non-linear analysis is also essential in reaching the conclu- 
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sion (c/. Figs. [5] and [3]) that, for Ji ^ J2, the pairing am- 
plitudes for several competing nodeless and nodal pairing 
channels are comparable. This last result is particularly 
important for experiments in the iron pnictides. Indeed, 
our result has anticipated the recent experimental find- 
ings |36[I37| that pnictide superconductors with nodeless 
or nodal gaps have a comparable maximum Tc. Our s 2^2 

X y 

is±), ^x^+y^ states appear to be consistent with ARPES 
measurements that find full gap at hole pockets. How- 
ever in contrast to s 2 ^2 , which is fully gapped on all the 
Fermi pockets as long as the pockets are not too large, 

Ai 

s^2^y2 state has nodes on the electron pockets. This may 
be the reason for seeing nodal behavior in P-doped 122 
compounds in contrast to the fully gapped behavior in 
K-doped 122 compounds. 

To summarize, we have shown that magnetic frustration 
effects lead to quasi-degeneracy among different pairing 
states. With the bandstructures appropriate for the iron 
pnictides, an extended Aig state is the dominant pairing 
state, but this state contains both Sx^y2 and Sx^^y2 com- 
ponents. Moreover, a Big state has a close-by ground 
state energy. The quasi-degeneracy makes it likely that 
the iron pnictides of different material families, or differ- 
ent dopings, have different superconducting states. Our 
detailed phase diagram contains a low-temperature phase 
with time-reversal-symmetry breaking, which also can be 
tested by future experiments. 
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